Finite range corrections near a Feshbach resonance and their role in the 

Efimov effect 



P. Dyke, S. E. Pollack, and R. G. Hulet, 

Department of Physics and Astronomy and Rice Quantum Institute, 
Rice University, Houston, Texas 77005, USA 
(Dated: Febmary 4, 2013) 

Abstract 

We have measured the binding energy of ^Li Feshbach molecules deep into the non-universal regime 
by associating free atoms in a Bose-Einstein condensate by modulating the magnetic field. We extract the 
scattering length from these measurements, correcting for non-universal short-range effects using several 
different methods. We find that field-dependent effective range corrections agree well with the data. With 
this more precise determination of the scattering length vs. field we reanalyze our previous data on the 
location of atom loss features produced by the Efimov effect Uj] and investigate effective range coixections 
to universal theory. 
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Efimov showed more than 30 years ago that three particles interacting via resonant two-body 
interactions could form an infinite series of three-body bound states if the two-body 5-wave scat- 
tering length a could be varied [|2|]. In the limit of zero-range interactions, the ratios of scattering 
lengths corresponding to the appearance of each bound state was predicted to be a universal con- 
stant, equal to approximately 22.7. The only definitive observations of the Efimov effect have been 
in ultracold atoms, where the ability to tune a via a Feshbach resonance [|3i |4|] has proven to be 
essential. Since the first evidence for Efimov trimers was obtained in ultracold Cs Isl, experiments 
have revealed both three- and four-body Efimov states in several atomic species. Although the Efi- 
mov effect has now been confirmed, several open questions remain, including a full understanding 
of the role of non-universal finite range effects. Accurate comparisons with theory require that 
these non-universal contributions be quantitatively determined and incorporated. 

We previously characterized the F = \,mp = \ Feshbach resonance in 'Li, which is located at 
approximately 737 G, by extracting a from the measured size of trapped Bose-Einstein conden- 
sates (BEC) assuming a mean-field Thomas-Fermi density distribution [lli|6|]. These data were fit 
to obtain a{B), the function giving a vs. magnetic field, which was used to assign values of a to 
Efimov features observed in the rate of inelastic three- and four-body loss of trapped atoms [UJ]. 
More recently, two groups have characterized the same Feshbach resonance by directly measuring 
the binding energy, E^,, of the weakly-bound dimers on the a > side of the Feshbach resonance 
|7-^. These measurements disagree with our previous measurements based on BEC size. The dis- 
agreement in the parameters characterizing the Feshbach parameters is sufficiently large to affect 
the comparison of the measured Efimov features with universal theory. 

In this paper, we report new measurements of Eb, which we fit to obtain a{B). The measurement 
of has fewer systematic uncertainties than the BEC size measurement, which is affected at large 
scattering length by beyond mean-field effects and by anharmonic contributions to the trapping 
potential. The extraction of a from can therefore be more accurate, and unlike the condensate 
size measurement, E^ is related to a for both thermal gases and condensates. We have measured E^ 
far enough from the Feshbach resonance that E\^ no longer depends quadratically on the detuning 
of B from resonance, as expected in the universal regime flSl |4|| . We show that a commonly adopted 
correction for non-universal finite range effects, which depends on a single value for the effective 
range, does not fit the data as well as field-dependent effective range corrections, or a more exact 



formalism arising from a two-channel model [llOll . We use the newly measured and corrected a{B) 



function to reanalyze our three- and four-body loss data to obtain improved locations of the Efimov 



features. 

Our experimental methods for producing EEC's and ultracold gases of ^Li have been described 
in detail previously [|6|]. Atoms in the \F = l,mf = 1) state are confined in an optical trap 
formed from a single focused laser beam with wavelength of 1 .06 jum. A bias magnetic field, 
directed along the trap axis, is used to tune a via the Feshbach resonance. For the new data 
presented here the axial and radial trapping frequencies are 4.7 Hz and 255 Hz respectively. We 
adjust the magnetic field to give a ~ 200 ao, where ao is the Bohr radius, and use forced evaporation 
to produce either ultracold thermal clouds with temperatures of ~1.5/iK, or condensates with 
condensate fraction that we estimate is greater than 85%. We then adiabatically ramp the field 



to the desired value and employ in situ imaging, either polarization phase-contrast [lllll when the 
density is high, or absorption imaging in less dense clouds. 

Atoms are associated into Feshbach molecules by resonantly oscillating the magnetic field at 
a frequency hv^od = £'b + £'kin, where Eh is taken to be positive for a bound state, and ^kin is the 



relative kinetic energy of the atom pair [|l2[|l3|l . The weakly -bound dimers formed in this way are 



lost from the trap through coUisional relaxation, presumably into deepl y-b ound vibrational levels 



12]. This technique has been used in studies of both homonuclear jl 



14 



15ll and heteronuclear 



lll6l - ll8n Feshbach resonances, in addition to the specific hyperfine state of 'Li studied here [|7|-|9|]. 

The oscillating field is produced by a set of auxiliary coils that are coaxial with the bias coils 
producing the Feshbach field. The amplitude of this field ranges from 0.1 G to 0.6 G and the 
duration of modulation ranges from 25 ms to 500 ms depending on magnetic field. The number 
of remaining atoms are measured as a function of the frequency of the oscillating field. We find 
that in the case of a BEC, the loss spectra are fit well by a Lorentzian lineshape. In the case of a 
thermal gas we fit the loss spectra to a Lorentzian convolved with a thermal Boltzmann distribution 

iBQQ]. 

Figure [H shows characteristic loss spectra at 734.55 G (where a ^ 1 100 ^o) for several different 
temperatures and modulation amplitudes. The Lorentzian component fits to a linewidth of 8 kHz, 
which provides a lower bound on the lifetime of the molecular state of 20yus. The temperature 
extracted from the thermal broadening of the loss feature of the 1 .5 yuK data driven with the smallest 
modulation amplitude of 0.14 G agrees within 20% of the measured temperature found by fitting 
the atomic distribution to a Gaussian. There is no observed thermal broadening of the BEC loss 
data, as it fits well to a pure Lorentzian. There is no systematic shift in the resonance location 
with temperature or modulation amplitude, but for large amplitude modulations and sufficiently 
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FIG. 1. (Color online) Magneto-association induced loss at B = 734.55 G, where a - llOOao. Main plot 
shows loss spectra for thermal gases with the following modulation amplitudes and temperatures: (a) 0.57 G 
and 20 //K, (■) 0.14G and 1.5 //K, and (•) 0.57 G and 1.5 //K. The solid curves are fits to Lorentzians 
convolved with themial Boltzmann distributions. The inset (♦) corresponds to a BEC with a modulation 
amplitude of 0.14 G. The solid black line is a Lorentzian fit to the condensate resonance and the vertical 
dashed line in the main figure is the resonance location E\^lh = 450 kHz found from this fit. 



low temperature we observe a nonlinear resonance at I 



similar nonlinear response was reported previously [llSll . 



E\,/h. No other subharmonics are seen. A 



The results of the measurement of binding energy vs. B are displayed in Fig. |2la). In the 
universal regime (see Fig. l2lb)), where a is much larger than any characteristic length scale of the 
interaction potential, = Jf/ma^, where m is the atomic mass [13i|4|]. The solid lines in Fig. |2] 
show the results of fitting in this universal regime to a, where a is given by the usual Feshbach 
resonance expression 

and where a\,g is the background scattering length, A is the width of the resonance, and 5oo is the 
location of the resonance. In the large a universal regime, the first term in Eq. (H) is negligible, and 
the fit has just two free parameters, and the product a^g A (values given in TableU]). In this case, 



E\j oc {B - Boof. In Fig.[2l;c) and (d), the same data is recast in terms of 7 = (mEb/ff-y^^, where 
for large a, y vs. 5 is a straight line, as is shown in Fig. |2td) . We find no detectable difference in 
£'b between a BEC or a thermal gas. 

Figure|2]suggests that the universal regime extends down to ~725 G, or ~ 1 2 G below resonance. 
Significant discrepancies between the measured E^, and universal theory are observed as the field 
is decreased further. It is not surprising that the universal regime is only a small fraction of A (174 
G), since the 'Li resonance is known to be intermediate between closed-channel and open-channel 
dominated, as the resonance strength parameter i'res - 0.56 I19I1 is neither «; 1 nor » 1 [4]. For a 
more precise determination of a it is desirable to extend the analysis into the non-universal regime, 
where short range attributes of the potential become appreciable. A simple two-channel approach 
to correct for finite range effects, suggested in Ref. [4] and applied to 'Li in Refs. [Q, 0], is to 
replace the universal binding energy expression with 



m 



{a-a+R*y' 



(2) 



where a = 31 ao is the mean scattering length n2(M (closely related to the van der Waals radius 



civdw = 32.5ao), and R* = al s^e^ = 55 uq is related to the resonance width [12 ill . The bound state has 
predominately open channel character only for a » AR* , which is the expected range of validity 
of Eq. |2] [|41]. The best fit to the data using Eqs. [T]and|2]is plotted in Figs. [2l 

Although Eq.|2]gives a somewhat better fit to the data than the universal binding energy relation, 
it clearly fails to represent the entire range of measurements. Higher order corrections to this 



theory offer little improvement to the overall fit quality [|4l |22||. A more complex solution to the 



two-channel model is given by Eq. 26 in Ref. [llOh . This solution is given in terms of R* and 



another parameter b, which is related to the van der Waals length and hence, to a. In order to agree 
with Eq. |2]in the proper limits, we identify b = ^a. Using the previously specified values of R* 
and a, the best fit to the data is shown in Figs.[2l The expected improvement over the simple model 
(Eq. O is borne out, as its range of validity extends to larger detunings from resonance. 

The simple two-channel approach represents the effective range of the potential, R^., with a 
single value. Since the 'Li resonance is not open-channel dominated, however, R^ exhibits con- 
siderable field-dependence over the width of the resonance. Properly accounting for this field 
variation should provide a better correction for finite range effects. In order to obtain Re(B) we 
numerically solved the full coupled-channels equations using realistic model potentials for both 
the singlet (closed channel) and triplet (open channel) potentials of the electronic ground state 
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FIG. 2. (Color online) Results of modulation spectroscopy using condensates (•) and 1.5 //K thermal clouds 
(■). (a) and (b): Ey^ vs. B, and (c) and (d): the same data plotted as y = {mE\,/H^y^'^ vs. B. The vertical 
eiTor bars corresponds to fitting uncertainty, while the horizontal error bars are the shot-to-shot variation of 
the magnetic field. The solid (black) curves are fits to Eq. ([1} using the universal expression E\, - H^/ma^; 
the dashed (green) curves are fits using the simple two-channel model, Eq. Q, with the parameters given in 
the text; the dot-dashed (blue) curves use the solution to the complex two-channel model given in Ref. lioh 
(Eq. 26), again with parameters as given in the text; and the dotted (red) curves incorporate the field- 
dependent effective range given in Fig.|3] using Eq. (UJ. The fits exclude data below 725 G, where a ^ 250 uq 
is no longer much greater than R*, and the validity of the non-universal corrections becomes questionable. 
The resulting Feshbach resonance parameters for the universal model and for the field-dependent effective 
range are given in Table |I] The fitted parameters for the latter are equal to those obtained from the two- 
channel models, to within experimental uncertainty, for fits constrained to B > 725 G. 
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FIG. 3. (Color online) Coupled channels calculation of a and for the F = l,mf = 1 Feshbach resonance 
in ^Li. The horizontal and vertical dashed (blue) lines indicate a = and Boo, respectively. /?e was fit to a 
polynomial expansion in the scaled field b - B - 12>1 .1 G to obtain Re(ao) = -6.2 + 3.50b - 9.2 x lO'^b^ - 
6.5 X lO-^Z;^ + 5.7 x 10''^ b\ 



of Li [1231. 12411 ■ These potentials have been refined by adjusting parameters, such as the potential 



depth and the shape of the inner wall, to give quantitative agreement with experimentally known 
quantities, which are primarily the locations of Feshbach resonances [6, 8], zero crossings and 



the binding energies of the least bound triplet molecule n25\\ . The scattering length and effective 
range are determined from the energy dependence of the 5- wave phase shift 



1 1 , 
kcotdoik) = — + -Rek + 
a 2 



(3) 



where frk~/m = E]^„. Figure |3] shows both a and R^. near the Feshbach resonance at 737 G. 
There is considerable variation in R^ over the width of the resonance, contrary to the assignment 



Re = -2R* nil 



or Re = 2{a - R*) [|4|]. In comparison, the coupled-channels calculation gives 
Re(Bco) - -6 Gq. Given this variation, it is not surprising that Eq. |2]does not describe the data well. 
The relation between the binding energy of a weakly bound state, or equivalently y, and a and 



Model 


Boo (G) 


A(G) 


fl^bg («()) 


dbg A (G ao) 


Universal 


737.67(10) 


-174 


-20.6 


3584(60) 


Coupled-channels /?e 


737.67(10) 


-174 


-19.8 


3450(60) 



TABLE I. Feshbach resonance parameters obtained by fitting Eb data to Eq. [T] using either the universal 
model (y oc 1/a) or using Eq. |4]with the coupled-channels calculation of the field-dependent values of R^. 
There are large uncertainties in Obg and A separately, but their product is well-defined by the data. The 
choice of A was guided by the coupled-channels calculation. The quoted uncertainties reflect shot-to-shot 
variations in the field and fitting uncertainties. Since the range of validity of the non-universal models is 
guaranteed only for ■« a, data below 725 G is excluded from the fit. In this field range, fits to the 
two-channel models, the simple one given by Eq. [2] l^, and the more complex solution given by Eq. 26 
of Ref. H, give parameter values that are the same, to within the uncertainties, as the field-dependent R^ 
model. The parameters from this model, given in bold, are our recommended values. 



i?e is given by y = l/a + ^Rej^ 



physical behavior for \Re/a\ «; 1 [|lOL 



27 



21, 



lis quadratic equation has one solution with the correct 



1- A 1- 



2R, 



(4) 



For \Re/a\ 1, y ^ l/a (l + ^Rg/aj. Figures |2c) and (d) show the results of fitting the measured 
values of y to Eqs. ([T]) and The agreement between theory and experiment is very good over 
a much larger range of the measurements, and we use this fit to define the Feshbach parameters, 
which are indicated in bold in Table H Since a » la^gl for all of the E^, data points, the data is 
insufficient to separately extract both abg and A. Given the precise knowledge of the location of 
the field where a = 0, Bo = 543.6(1), found in our previous work [|6|], a logical choice would be 
to fix A = 5o - Boo = -194.1 G. The coupled-channels calculations, however, offer a useful guide 
to assigning values to these parameters separately. We find that A = -174G gives a somewhat 
better fit to the coupled channel results over a larger range in a, so we adopt this value. We stress, 
however, that the fit to the data strongly constrains the product abg A, but not each parameter 
separately. A similar procedure was followed in Refs. \SAm- Our Feshbach parameters agree with 
Ref. where they find Boo = 737.8(2) G, but slightly disagrees with Ref. [js!], which reports 
5oo = 738.2(4) G. 

We now turn our attention to the three- and four-body Efimov features previously reported in 
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Ref. yj]. Figure m shows the measured three-body loss rate coefficient L3 plotted vs. a, where 
the correspondence between measured values of B is now determined by the new Feshbach pa- 
rameters given in Table H While L3 generally scales as a"* it is punctuated by several minima and 
maxima, which are due to Efimov molecular states. The previously reported Efimov maximum 
a^, corresponding to the second Efimov trimer, is no longer visible because the upward shift in 
the resonance position by 0.7 G relative to the previous measurement [ li places this feature in the 



regime where the loss rates are limited by quantum mechanical unitarity [|29L 13011 . For the same 
reason, the effect of the second Efimov tetramer associated with the second trimer 2) is also not 
visible in measurements of the four-body loss rate coefficient L4 (not shown). The fitted locations 
of the remaining features are given in the second column of Table HH where the caption provides a 
key to the notation. Not all of the features given in Table|II]are indicated in Figure HI but expanded 
views of both L3 and L4 are found in Ref. i^. 



The origins of three of the features in Table III indicated by square brackets, are uncertain. 
The feature a*2 is nominally located at the atom-dimer resonance where the energy of the second 
Efimov trimer merges with the atom-dimer continuum. Relatively sharp peaks in L3, located near 

n hn 

the expected atom-dimer resonance, were previously reported for ^^K 113 8ll . and 'Li [llj,|39|]. Since 
a large dimer fraction is unexpected, a model was developed to explain the presence of enhanced 
loss even without a large population of dimers ^38^. In this model, each dimer produced in a 
three-body recombination collision shares its binding energy with multiple atoms as it leaves the 
trap volume due to the enhanced atom-dimer cross section [13 811 . Recent Monte Carlo calculations, 
however, conclude that the resulting peak from this avalanche mechanism is too broad and shifted 
to explain the observations [|40l1 . The remaining two features, a*2 ^ and nominally located 

at dimer-dimer resonances, where the energy of a tetramer merges with the dimer-dimer threshold 
ll36ll . Their assignment also remains tentative, since their observation requires a significant and 
unsubstantiated dimer population. 

The third column in Table pi g ives the predictions of universal scaling. Pioneering scaling 
predictions for the three-body [31| and four-body sectors [3^ 371, have largely been replaced by 
the more precise theoretical determinations cited in Table UIl Four significant digits are given to 
reflect the stated precision of the scaling relations. Since the relative positions of all features are 
expected to be universally connected, the position of only one is needed to completely fix the 
remaining. We choose the recombination minimum of the second trimer, , for this purpose as 





Feature 


Experiment (ao) 


Universal Scaling in a (ao) 


Universal Scaling in y (^o) 


a> 


a| 


88(4) 


61.78* 


17.57 




a- 


1402(100) 


(1402) 


(1402) 




14 


426(20) 


313.5*" 


293.8 




\a* 1 

[2,1 J 


919(50) 


688.5^ 


680.4 






1902(200) 


2127^ 


2132 


a <0 




-241(8) 


-294.3+ 


-277.7 




jL 




-6679^^ 


-6716 




T 


-93(4) 


-125.2** 


-75.40 




T 
<2 


-235(10) 


-268.5** 


-249.6 




T 
^2,1 


-4040(700) 


-2841** 


-2857 




T 
^2,2 




-6094** 


-6129 



References: * 1311]; + ISJ; * 131; ^ ISJ; ^ 1351 



TABLE II. Locations of Efimov features, given in units of ao, of the three- (L3) and four-body (L4) loss 
coefficients. The experimental values of a are extracted from the measured fields using Eq. {[T}. The hor- 
izontal lines indicate features that are too near the resonance to be observed. The estimated uncertainties 
include fitting uncertainties, as well as the uncertainties in the Feshbach parameters. For a > 0, at de- 
notes the recombination minimum of the i* Efimov trimer. The origins of three features, labeled as |^a* j, 
^a*2 j], and [^a* 2], have not been identified, but roughly correspond to expected locations of atom-dimer and 
dimer-dimer resonances. For a < 0, denotes the Efimov resonance where the i* trimer merges with the 
free atom continuum. The remaining features, aj j, arise where the j"' tetramer associated with the i"^ trimer 
merges with the free atom continuum. The final two columns give the predicted locations of the features 
using universal scaling in a and y, respectively. The universal scaling relations were obtained from the 
indicated references. The scaling is anchored by the measured location of a^, which, as an input, is denoted 
by parentheses. For the last column, the scaling relations are applied to ratios of y rather than a, and the 
table entries are the values of a corresponding to these y. 



10 



1 1 1 1 I I I — 1— 




10^ 10^ 10"^ 







10"^ -10^ -10^ 



FIG. 4. (Color online) Three-body loss rate coefficient L3 vs. a. The points are values extracted from the 
measured trap loss, with the green, blue, red, etc correspond to different trap frequencies, as reported in 
Ref. Uj]. The sohd line shows universal scahng in a, where the positions of the features are determined by 
the single feature a}^. The dashed lines show the a' scaling in the absence of the Efimov effect. The fitted 
width parameters are 77+ - 0.038 and ?/_ =0.12, for the t? > and a <0 sides of the resonance, respectively. 



it is a well-defined feature that occurs at sufficiently large a (~1400(3o) to be insensitive to short- 
range effects, while also being small enough in magnitude to not be hypersensitive to B. While the 
agreement with universal theory of ~30% is good overall, some of the features, in particular a[ and 
the lowest tetramer a[j, occur in the non-universal regime where \Rela\ is not small. We attempted 
to correct the universal theory for the effect of finite range using the same strategy applied to the 
dimer binding energy, that is by applying universal scaling in y"^ (Eq. |4l) rather than a. To lowest 
order, the correction to 1 la is ^R^-fa^. The fourth column in Table |ll] gives these values. We find 
that such a replacement improves the agreement with experiment for features on the a < side of 
the resonance, but for a > the agreement is actually made worse. For the a > side of resonance, 
we note that a correction of opposite sign, corresponding to replacing with |i?el, gives the best 

11 



agreement. For example, the feature a| observed at a = 88 ao is predicted to be located at 83 ao 
under this substitution. The fact that this sign change improves the fit is an empirical observation 
for which we have no theoretical explanation. We note that an eff'ective field theory for short range 
interactions has been developed in which corrections to universal scaling of three-body quantities 
are computed to 0(Re) and that they have been applied to the F = l,mf = Feshbach resonance 
in 'Li [14 in . It would be interesting to apply the same analysis to the F = l,mf = 1 resonance to 
compare with the data presented here. 

A measure of universality across the Feshbach resonance may be obtained by evaluating the 
ratio a^la^. Universal scaling in a implies a^la'^^ = -4.76 1321], whereas experimentally, we 
find a^l a\ = -5.82. We disagree with a previous measurement for the \F = l,mp = 1) state in 
'Li, which found d^l a\ = -4.61 |18i|420. The discrepancy can be attributed to difi'erences in the 



measured Feshbach parameters described previously. In addition, the Efimov features observed in 
Refs. [V.Is] are not as sharp as those reported here, and this may affect the precision for which the 
location of a feature is extracted. The width is quantified by the fit parameter 77, which is related 
to the lifetime of the Efimov molecule pll] . For Ref. [j8|], 77+ = 0.17 and rj^ = 0.25, corresponding 
to the a > and a < sides of resonance, respectively, while we find 77+ = 0.038 and //_ = 0.12. 
These large differences, at least in case of 77+, may indicate that 77, and hence the dimer lifetime, 
has an interesting and unexpected temperature dependence, since the a > data in Refs. 13, ISQ 
is obtained with a thermal gas, while in our experiment the gas is cooled to nearly a pure Bose 
condensate. 

It was pointed out recently that the location of the first Efimov trimer resonance when scaled 
by the van der Waals radius a^dw, is remarkably similar for multiple unconnected resonances in 
the same atom [|43|] . as well as for different atomic species [|44|l. These observations suggest that 
there is no need for an additional "three-body parameter" to pin down the absolute positions of the 
Efimov features, but rather, that this scale is set by short-range two-body physics [|44l-l48l]. For the 
measurements reported here, -a7 /a,wvv = 7.42(3), which is reasonably close to the range of 8-10 



reported in most other cases [|43 



4411 



The determination of the Feshbach parameters for 'Li in the \F = \,mp = 1) state by di- 
rect measurement of the dimer binding energy is a significant improvement over our previous 
measurement using condensate size. Using these more precise parameters we find that the over- 
all agreement between the experimentally determined locations of three- and four-body Efimov 
features and those obtained from universal scaling is in the range of 20-30%. Since we use the 



12 



location of only one feature as input, the agreement strongly supports the contention of universal 
scaling across the Feshbach resonance. The origin of features nominally located at atom-dimer 
and dimer-dimer resonances remains an open question. A direct measurement of the equilibrium 
dimer fraction could help to resolve this issue. 
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